arXiv: 1505.06714v2 [nucl-th] 8 Feb 2016 


EPJ manuscript No. 

(will be inserted by the editor) 


Quark matter under strong magnetic fields* 

Debora Peres Menezes^ and Luiz Laercio Lopes^’^ 

^ Depto de Ffsica - CFM - Universidade Federal de Santa Catarina - Florianopolis - SC - CEP 88.040 - 900 - Brazil 
^ Centro Federal de Educagao Tecnologica de Minas Gerais, Campus VIII - Varginha - MC - CEP 37.022-056 - Brazil 

Received: date / Revised version: date 

Abstract. We revisit three of the mathematical formalisms used to describe magnetized quark matter in 
compact objects within the MIT and the Nambu-Jona-Lasinio models and then compare their results. The 
tree formalisms are based on 1) isotropic equations of state, 2) anisotropic equations of state with different 
parallel and perpendicular pressures and 3) the assumption of a chaotic field approximation that results in 
a truly isotropic equation of state. We have seen that the magnetization obtained with both models is very 
different: while the MIT model produces well-behaved curves that are always positive for large magnetic 
helds, the NJL model yields a magnetization with lots of spikes and negative values. This fact has strong 
consequences on the results based on the existence of anisotropic equations of state. We have also seen 
that, while the isotropic formalism results in maximum stellar masses that increase considerably when the 
magnetic fields increase, maximum masses obtained with the chaotic field approximation never vary more 
than 5.5%. The effect of the magnetic field on the radii is opposed in the MIT and NJL models: with both 
formalisms, isotropic and chaotic field approximation, for a Hxed mass, the radii increase with the increase 
of the magnetic field in the MIT bag model and decrease in the NJL, the radii of quark stars described by 
the NJL model being smaller than the ones described by the MIT model. 

PACS. 24.10.Jv Relativistic models - 21.65.Qr Quark matter 


1 Introduction 

The influence of magnetic fields (B) in the formation, con¬ 
stitution and evolution of stars is a topic of intense inves¬ 
tigation and discussion. Some astronomers believe that 
magnetism plays an essential role in the star formation, 
second only to gravitation [T]. As the star burns its fuel, 
depending on its mass, it may end as a white dwarf or 
a neutron star and in both cases, magnetic fields can be 
quite intense. Neutron stars are very dense compact ob¬ 
jects with magnetic fields generally of the order of 10^^ G 
and they are mainly detected as pulsars, powered by ro¬ 
tation energy, or as accreting X-ray binaries, powered by 
gravitational energy. However, since 1979, other classes of 
neutron stars have been observed as either soft gamma-ray 
repeaters or as anomalous X-ray pulsars with no binary 
companion. These objects bear magnetic fields of the or¬ 
der of 10^"^ — 10^® G and were named magnetars [2]. 

If one observes the QGD phase diagram, neutron stars 
lie in a region of low temperatures and high baryonic 
chemical potentials fi (or baryonic densities, if one prefers). 
If these objects are affected by the presence of strong 
magnetic fields, other regions of the QGD phase diagram 
should be affected as well 0, i. 

The region of high temperatures and low chemical po¬ 
tentials is the relevant regime for the study of heavy ion 
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collisions [S] and in these experiments, the field intensity 
probably decreases very rapidly, lasting for about 1-2 fm/c 
only, but a possible signature of the strong magnetic fields 
has been investigated 0, H, 0,0, m- Moreover, this 
region of high-T/low-/i, has already been exploited using 
Lattice (LQCD) simulations. The QGD phase diagram 
presents a first order transition at low-T/high-/i region 
and a crossover region at high-T /low-fj, and these two re¬ 
gions may be connected by a critical end point SI, m- 
LQCD results point out, in accordance with most model 
predictions, that the crossover observed at H = 0 does not 
disappear when B ^om,m, m, sa- Nevertheless 
the region of low-T/high-/i remains unavailable to LQCD 
assessments and hence, the physical properties related to 
this region can only be investigated with the help of (ef¬ 
fective) models. 

Other important aspects related to magnetized quark 
matter are the richness of its internal substructures and 
the inverse magnetic catalysis (IMG). Although the results 
are model dependent, the internal structure of the QGD 
phase diagram has been shown to be very rich m , pre¬ 
senting lots of intermediate phases, due to the small jumps 
in the quark dressed masses, which are related to the filling 
of the Landau levels. Concerning the IMG, LQCD calcula¬ 
tions predict an inverse catalysis, with the pseudo-critical 
temperature decreasing with B [niin], while most effec¬ 
tive models predict an increase of the pseudo-critical tern- 
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perature with B HZ], m- It is a common belief that this 
deficiency in effective models is due to the fact that they 
do not account for back reaction effects (the indirect in¬ 
teraction of gluons with the magnetic field) [Ifij . 

One important point that should not be disregarded 
refers to the contribution of the electromagnetic interac¬ 
tion to the pressure(s) and energy density, a term pro¬ 
portional do B^, where B is the magnetic field strength. 
Astrophysicists advocate that the magnetic field in the 
surface of magnetars is of the order of fO^^G. However, 
according to the Virial theorem one could expect fields 
three or four orders of magnitude stronger in their inte¬ 
rior. To cope with this difference in strength, an ad hoc ex¬ 
ponential density-dependent magnetic field was proposed 
in ref. m and adopted in many other works ED, ED, 

m ED, ED, m, ED, ED: ESI: EQ], m. M- Another 

similar prescription for an energy density dependent mag¬ 
netic field was proposed in Ref. ESj • The second prescrip¬ 
tion seems more natural because the quantity used in the 
TOV equations EH to calculate the macroscopic quan¬ 
tities is the energy density and not the number density. 
Moreover, in E2 ]j it is shown that this choice reduces the 
number of free parameters and it is given by: 

H(e) = Ho(^y+ (1) 

where Bq is the fixed value of the magnetic field, cq is the 
energy density at the center of the maximum mass neutron 
star with zero magnetic field, cm is the energy density 
of the matter alone, avoiding a self-generated magnetic 
field, and 7 is any positive number. In the present work 
we take 7 = 8 and thus, we reduce the number of free 
parameters from two, used in the prescription proposed 
in m to only one. Nevertheless, as pointed in ref. ESj: 
for chaotic magnetic fields, Eq. o yields a parameter-free 
model, while for the standard density-dependent magnetic 
field, the macroscopic properties are strongly dependent 
of these two non-observables parameters. Another point 
worth mentioning is that the magnetic field is no longer 
fixed for all neutron star configurations. Each equation of 
state (EOS) produces a different value for cq that enters 
in Eq. (HD. In this sense, eg is also another parameter, but 
it is not arbitrary as 7 and the two parameters that can 
be arbitrarily changed in the parametrization proposed in 
ref. EO] • It comes directly from the model that defines the 
EOS, as the baryonic density, pressure and energy density 
that enter as input to the TOV [51] . 

One should bear in mind that these two ansdtze, vi¬ 
olate Maxwell equations because the divergent of B is 
no longer zero. To minimize this problem, in the present 
work only the contribution to the term proportional to 
B^ will be taken as density dependent. The EOS obtained 
from magnetized matter will always carry a fixed B value. 
With this choice, if we apply the models to describe stellar 
matter, we guarantee thermodynamical consistency in the 
EOS and still force matter at the surface of the star to be 
subject to a magnetic field that is not too strong. It is also 
important to stress that if a density dependent magnetic 
field is used through out the calculations, numerical re¬ 


sults are almost coincident with the ones we present next 
but, besides violating Maxwell equations, thermodynami¬ 
cal consistency would also require an extra term, normally 
disregarded in the literature. This point will be specifically 
emphasized when the equations are displayed in the next 
sections. 

In the present work we study quark matter possibly 
present in strange (quark) stars. According to the Bodmer- 
Witten conjecture ES, El]> EZ]: the interior of a neutron¬ 
like star does not consist primarily of hadrons, but rather 
of quark matter composed of deconfined up, down and 
strange quarks, plus the leptons necessary to ensure charge 
neutrality and /3-equilibrium. When a model is chosen, it 
is important to verify if it satisfies the necessary condi¬ 
tions that ensure that quark matter is the true ground 
state: two-flavor quark matter must be unstable (i.e., at 
zero temperature its energy per baryon has to be larger 
than 930 MeV, the iron binding energy) and the three- 
flavor quark matter must be stable (i.e., its energy per 
baryon must be lower than 930 MeV, also at T = 0). 
Not all quark matter models satisfy these conditions and 
the NJL model only provides absolute stable matter when 
it is magnetized ES]- As already said, according to the 
Virial theorem and assuming uniform field and mass den¬ 
sity, the maximum magnetic field allowed in a gravitation¬ 
ally bound star is smaller than 10^® G. However, a quark 
star is self-bound by the nuclear interaction and a simple 
estimative of the maximum allowed magnetic field gives 
fields of the order of 10 ^° G ESIIIS]- We use these values 
as a guide to our study. 

It is worth mentioning that calculations that take into 
account the effects of the anomalous magnetic moments 
(AMM) of the quarks forsee that the maximum possible 
magnetic field strength is 8.6 xlO^^ G due to the fact 
that the u quark polarization would become complex for 
higher values m- On the other hand, works on hadronic 
stellar matter show that the influence of the (nucleonic 
and hyperonic) AMM on the EOS is minor El]) [H]) if the 
magnetic fields are of the order oi B = 10^® G. However, 
if H > 10^®G, the inclusion of the anomalous magnetic 
moment contribution stiffens the stellar matter EOS, and 
may originate a total spin polarization of neutrons [43] . In 
the present work, we do not consider the coupling between 
the field and the quark anomalous magnetic moment, but 
we think this problem should be tackled in a future work. 
We come back to this point in the conclusions of the paper. 

Despite the fact that usually a mean field approxima¬ 
tion is necessary when one uses effective models, restrict¬ 
ing their scope and the interpretation of their results, they 
can be very helpful if one wants to understand the physics 
not assessed by LQGD. Thus, different models have been 
recently used to describe quark matter subject to mag¬ 
netic field, the most common ones being the MIT bag 
model [U] and the Nambu-Jona-Lasinio model [IS]. In 
the present work we revisit both of them within three dif¬ 
ferent formalisms employed in the literature to describe 
stellar matter. 

The first formalism considers that quark matter is iso¬ 
tropic m and it was used in the seminal works on mag- 
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netars described either by a quark matter m or by a 
hadronic matter equation of state [48] and in subsequent 
calculations [32] , [35] . The first calculations involving quark 
matter described by the NJL model followed the same line 
of understanding [53], [33], [53], [53] , 

The second formalism is based on the fact that the 
EOS cannot be truly isotropic, once the components of 
the energy-momentum tensor are not equal, giving differ¬ 
ent contributions for the parallel and perpendicular pres¬ 
sure [33], [33], EH, [55], [S3]. Moreover, under strong mag¬ 
netic fields, the 0(3) rotational symmetry breaks, and this 
is another argument that supports the existence of pres¬ 
sure anisotropy. Once the EOS is obtained, one observes 
at what magnetic field the pressures start to deviate from 
each other and this is the limit usually taken as input to 
the Tolman-Oppenheimer-Volkoff equations (TOV) [33] . 
As pointed in ref. m, this limit is around 3.1 x 10^®G. 
In practice, the EOS used to compute the stellar macro¬ 
scopic quantities is isotropic [31], [55], [53] ■ However, this 
approach means that one recognizes that the system is 
anisotropic but ends up ignoring this fact. At least two 
recent works face this problem: in one of them |31] , the au¬ 
thors treated the anisotopic pressures as a perturbation in 
a way similar to the Hartle-Thorn method, generally used 
for slowly rotating neutron stars. In the other one [53] , an 
axisymmetric geometry is assumed and the Einstein equa¬ 
tions are solved with the adoption of a cylindrical symmet¬ 
ric metric. Older works also tackle this subject [55], [55]. 
In the present work we do not intend to reproduce these 
more sophisticated treatments, but we will investigate the 
effects of strong magnetic fields on the anisotopic pres¬ 
sures. 

Einally, the underlying assumption of the third for¬ 
malism is that in the presence of anisotropies, the concept 
of pressure has to be taken with care [57] [53]. Based on 
the concepts discussed in these two books, a small-scale 
chaotic field is used and the stress tensor is modified ac¬ 
cordingly, so that the resulting EOS is a truly isotropic 
one [33] . 

Next, a comparison of the results obtained with the 
three formalisms is shown and discussed. Although parts 
of this material are already available in the literature, a 
correct comparison can only be made if the same choice 
for the magnetic field is used. In the present work we use 
the energy density dependent magnetic field given in Eq. 
m to compute our results with the three formalisms. We 
would like to emphasize that our aim in the present work 
is not to justify, criticize or choose any of the formalisms 
just mentioned. We restrict ourselves to the analysis and 
comparison of the results. 

The present work is organized as follows: in Section (5] 
we present the EOS obtained for the MIT bag model in 
the presence of a magnetic field in one preferential direc¬ 
tion with the three possible formalisms just discussed: an 
isotropic EOS, an anisotropic EOS with different parallel 
and perpendicular pressures and another isotropic EOS re¬ 
sulting from the introduction of a chaotic field. The main 
results are then shown and discussed. In Section [3] the 
same steps are taken for the NJL model. In Section |31 


some of the results obtained with both models are com¬ 
pared. Einally, in the last Section, the final conclusions are 
drawn. 


2 The MIT Bag model 

The EOS for magnetized quark matter described by the 
MIT bag model has already been extensively studied [39] , 
[31], [51], [55], [53]. We next show only the Lagrangian 
density and the resulting EOS. We depart from the fol¬ 
lowing Lagrangian density: 

C = Cf+Ci-\F^.F>^\ ( 2 ) 

which contains a quark sector, Cf, a leptonic sector, Ci, 
and the electromagnetic contribution = d^A^ — dvA^. 
We use a static and constant magnetic field parallel to 
the z direction and hence, we choose the gauge = 
6 fj, 2 XiB and the energy levels in the x and y directions are 
quantized. 

The leptonic sector is described by 

= in [ 7 ^ (id^ - qiA^) - mi] ijji , (3) 

where I = e,fi. As we restrict ourselves to the T = 0 
case, the star has no more trapped neutrinos, which have 
escaped and carried a huge amount of energy while the 
star cooled down. 

The thermodynamical potential for the three flavor 
quark sector, 12 /, can be written as 

% = -Pf = - X! 

/ 

where Pf represents the pressure, £f the energy density, 
/i/ the chemical potential, and Uf the quark number den¬ 
sity. A similar expression can be written for the leptonic 
sector. 


2.1 MIT - Isotropic EOS 

In order to obtain the complete EOS, the pressure, and 
baryonic density were calculated in [35], [33], m and are 
given by 


Pf= E 

A;=0 


o^k- 


\Qf\BNe 




—s/(fc, B)^ In 



- Sfik,B)^ 
Sf(k,B) 


-s, 


(5) 


where Nc is the number of colors, g/ is the electric charge 
of each quark, B is the magnetic field strength and the 
quark masses are mu,d = 5 MeV, ms=120 MeV, B is the 
bag constant fixed as 148 MeV^/'^, oq = 1, cek>o = 2, 


Sf(k,B) = ^m'j + 2\qf\Bk 


( 6 ) 
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and 


with 


n = 


3 ’ 






E 


O^k 


27r2 


kpj, 


(7) 

( 8 ) 


where kpj = yM/ ~ Sf{k,B)‘^. At T — 0, the upper Lan¬ 
dau level (or the nearest integer) is defined by 


- mj k\j 
2\qf\B ^ 2\qf\B 


(9) 


and the energy density can be easily obtained from Eqs. 
(gl),® and ([5|). 

In the description of compact stars, both charge neu¬ 
trality and chemical equilibrium conditions have to be im¬ 
posed [59j . The first condition can be written for quarks 
and leptons as 


2n„ = nd + ns+i{ne + n ^), (10) 


and the second condition can be written as 


/^s — — fdu “f Me: Me — M/^: (71) 

where the lepton densities can be calculated through Eq. (|S]), 
with appropriate substitutions for the masses and elec¬ 
tric charge. The lepton masses are nie = 0.511 MeV and 
= 105.66 MeV. Eor the electron and muon pressure 
{Pi) and energy density (e/) we use Eqs. (O and the lep- 
tonic analogue of Eq. O, respectively, all with S = 0. In 
all cases, Nc = 1 for the leptons. 

The final expressions for the total pressure and energy 
density are given by 


Pf + Pi + B{e)y2, 

(12) 

ef + ei+B{e)^/2. 

(13) 


2.2 MIT - Anisotropic EOS 

The magnetization of the system is given by 
M = dP/dB, 


(14) 


and for the quark sector of the MIT bag model, it reads: 


^f,max 


E 


oik 


k=0 


kfWc 

dTT^ 


M/VM/ -s/(A:,B)2 


-{m'j + A\qf\Bk) In 



y^M/ - Sf{k,B)^ 
Sf{k,B) 


(15) 


In an anisotropic system, the parallel and the perpen¬ 
dicular components of the pressure can be written in terms 
of the magnetization, as [35]: gl]: |S1]: [33] : 


P|| = Pf + Pi-B{e)y2, Pp= Pf + Pi-MB + B{e)y2. 

(16) 

For a magnetic field in the z direction, the stress tensor 
has the form: diag(B^/2, i?^/2, —5^/2) and this explains 
the difference in sign appearing in the parallel and per¬ 
pendicular pressures. The energy density is the same as in 
the isotropic system, given by Eq. m- 

The magnetization for the lepton sector can be read 
off Eg. (11511 with the appropriate substitutions specified 
below. 


2.3 MIT - Chaotic field 

According to m, as the components of the stress tensor 
are not equal (they differ in sign), the quantity ±i3^/2 
cannot be simply added to the pressure terms. As already 
stated above, for a magnetic field in the z direction, the 
stress tensor is given by diag(i?^/2, 5^/2,—5^/2). One 
way of circumventing this problem is proposed in |58j . 
where it is shown that, for a small scale chaotic field, 
the concept of isotropic pressure is recovered because the 
stress tensor is diag(i?^/6, B^/6, i3^/6), the 0(3) rota¬ 
tional symmetry remains valid and hence, the pressure 
becomes: 


P = i < 7^* >= 1 
3 * 3 


p 2 



E 

6 ■ 


(17) 


Note that Eq. 113 represents the true thermodynamic 
pressure, since the components of the stress tensor are 
equal [SS], in opposition to the two first formalisms. We 
would like to point out that due to strong heat transporta¬ 
tion, the existing turbulence could create a chaotic field. 
Of course, as the star cools down, the fields could become 
ordered. The possible mechanisms responsible for creating 
and maintaining a chaotic field would have to be better 
investigated and are beyond the scope of the present work. 

Within this formalism, the expressions for the total 
pressure and energy density are given by: 

Pchao = Pf + Pi+ B{e)y6, (18) 

^chao = £/ + €[+ i?(e)^/2. (19) 


2.4 Results MIT 

We next show our results for the EOS obtained from the 
three formalisms just discussed. We have used cq = 6.93 
fm“^ in Eq. o because this value is the central energy 
density of the maximum mass obtained with the MIT 
model for non-magnetized matter with our choice of val¬ 
ues for the quark masses and the Bag constant. Before 
we proceed, it is important to stress that in Eqs. dm), 
m and m, the magnetic field is taken as constant in 
Pf and Pi and Eq.(|T|) is used in the last terms only. The 
same holds for the energy density terms. 
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Fig. 1. Magnetization as a function of the number density for 
the MIT bag model. 


Fig. 2. MIT bag model EOS obtained for B = 10*^® G without 
the term proportional to 


We first analyze how the magnetization varies with 
density for different values of the magnetic field in Fig. [TJ 
It is clear that the number of van alphen oscillations are 
larger for lower values of the magnetic field as expected, 
due to the larger number of filled Landau levels. As the 
magnetic field increases, the magnetization of the system 
also increases and hence, a stronger effect is felt on the 
anisotropy of the system, as can be seen in Eq. (HU). 

For B — 10^^ G, the magnetization is negligible and no 
visible difference can be seen when the perpendicular and 
parallel pressures are plotted. Moreover, the contribution 
coming from the term proportional do is also too small 
and, either it is included or not, the plots are coincident. 
For these reasons, we just plot results for magnetic fields 
larger than 10^® G. In Figs. (ED and m , we plot the parallel 
and perpendicular pressures for B = 10^® G respectively 
without and with the contribution proportional to the B^ 
term. In Figs. ED and (ED, O and ED, the same is dis¬ 
played for a larger magnetic fields, equal to B = 3 x 10^® G 
and B = 10^® G. By analyzing these figures, it is evident 
that the contribution of the B^ term has quite drastic ef¬ 
fects on the point where the pressures start to split and 
hence, as this term is important in stellar matter EOS, 
it is indeed difficult to justify the use of the TOV equa¬ 
tions for magnetic fields larger than 10^® G. One can also 
see, from Figs. ED, ED and ED that the isotropic pressure 
does not deviate much from the perpendicular pressure, 
the only difference being due to the magnetization of the 
system, which is still relatively small for these fields. The 
parallel pressure goes to zero at energy densities typical 
to neutron star core ii B = 10^® G, but for lower magnetic 
fields, the decrease starts at much higher densities. 

In Fig. (ED, the EOS obtained with the assumption of 
chaotic fields are shown. They deviate very little from non- 
magnetized matter, with a consequent small variation in 
the stellar properties that are computed next. 

We now use the EOS obtained from the three differ¬ 
ent formalisms we have discussed to compute the macro¬ 
scopic properties with the help of the TOV equations EU • 
In Figs. El uni ED we display the mass-radius relation for 



Fig. 3. MIT bag model EOS obtained for B = 3 x 10*^® G 
without the term proportional to B^ 



Fig. 4. MIT bag model EOS obtained for B = 10*^® G without 
the term proportional to B® 
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Fig. 5. MIT bag model EOS obtained for B = 10*^® G with 
the term proportional to 



Fig. 6. MIT bag model EOS obtained for B = 3x 10^® G with 
the term proportional to B^ 



Fig. 7 . MIT bag model EOS obtained for B = 10*^® G with 
the term proportional to B^ 



Fig. 8. MIT bag model EOS for different values of the chaotic 
magnetic field. 



Fig. 9. MIT bag model: Mass-radius relation within isotropic 
pressure for different values of magnetic field. 



Fig. 10. MIT bag model: Mass-radius relation within 
anisotropic pressure for different values of magnetic field. 
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Fig. 11. MIT bag modeliMass-radius relation within chaotic 
magnetic field approximation for different values of magnetic 
field. 


Table 1. Properties of the maximum mass quark star 


Model 

B 

■^max 


R 



(G) 

(Me) 

{Mq) 

(km) 

(fm-^) 

M^ITisotropic 

TcF 

1.81 

2.32 

9.99 

6.77 

anisotropic 

TcF 

1.81 

2.32 

9.99 

6.83 

MITchaotic 

TcF 

1.81 

2.32 

9.99 

6.80 


M^ITisotropic 

10 "“ 

1.82 

2.33 

9.97 

6.77 

M^ITanisotropic 

TcP 

1.82 

2.33 

9.95 

6.85 

MITchaotic 

TcP 

1.82 

2.32 

10.03 

6.63 


M^ITisotropic 

3 . 10 "“ 

1.89 

2.42 

9.78 

7.52 

M^ITanisotropic 

3 . 10 "“ 

1.85 

2.36 

9.58 

7.63 

MITchaotic 

3 . 10 "“ 

1.83 

2.33 

10.08 

6.63 


M^ITisotropic 

10 "“ 

2.16 

2.90 

10.68 

6.33 

M^ITanisotropic 

TcP 

1.99 

2.49 

9.31 

8.39 

MITchaotic 

TcP 

1.99 

2.66 

10.97 

5.87 


isotropic, anisotropic and chaotic magnetic field approxi¬ 
mation respectively, and in Table [T] we display the prop¬ 
erties for the maximum mass quark star. Had we chosen a 
smaller value for the strange quark mass, we would have 
obtained larger maximum masses, as can be seen, for in¬ 
stance, from Fig. 5b in Ref. [SD]. It means that a 2 Mq 
star could be attained with another choice of parameters, 
but our general conclusions remain the same. The same 
statement is valid with respect to the choice of the bag 
constant B fixed as 148 MeV^/'^, a value that satisfies the 
Bodmer-Witten conjecture once the stability window is 
investigated and gives a maximum stellar mass not too 
low, as can be seen also from Fig. 5b in Ref. m- 

The maximum mass only shows a real increase as com¬ 
pared with the one obtained from non-magnetized matter 


(which is identical to the results obtained for B = 10^^ G) 
for H > 3 X 10^® G for an isotropic EOS and for B ~ 10^® 
G if the chaotic field or the anisotropic pressure is used. 
This is easily understood because the EOS for chaotic field 
approximation is softer than the one for the isotropic EOS, 
as one sees by comparing Eqs. m and (fTSl) . In solving 
the TOV equations, we have used the perpendicular pres¬ 
sure (and not the parallel one) given in Ea. dTBl) because 
the parallel pressure goes to zero at certain energy densi¬ 
ties, as discussed above. Hence, for the magnetic fields of 
interest, the maximum mass would be too low. This effect 
was also shown in m, where the authors have computed 
the TOV equations with both pressures. Of course, it is 
very difficult to give a reasonable physical interpretation 
for this kind of calculation that only considers one of the 
existing pressures, but we present it here for the sake of 
completeness. 

Looking at the results shown in Table [U one can see 
that the maximum stellar mass increases by approximately 
20% as compared with the non-magnetized star (equal to 
10^^ G if the isotropic EOS is used and only by 5.5% 
with the use of the chaotic field formalism. Our results 
come directy from the fact that the EOS obtained with 
the isotropic formalism is much harder than the one ob¬ 
tained with the chaotic approximation (it has an extra 
1/3 factor in the pressure of the magnetic field) for strong 
magnetic fields. The results obtained with the perpendicu¬ 
lar pressure of the anisotropic EOS are similar to the ones 
obtained with the chaotic field approximation however, al¬ 
ways with a smaller radii. Indeed we see that higher the 
magnetization, the more compact the quark star is. 

Another subject that we would like to discuss is related 
to the recent improvements of both, theory and observa¬ 
tions of neutron star radii. Based on chiral effective theory, 
the radii of the canonical 1.4M0 neutron star radius was 
constrained to 9.7 - 13.9 Km in [5T]. In [52], [53] a limit of 
12 km for the 1.4Mq was predicted, while in m this limit 
was set to 13.1 km. On the other hand, in [^ and |66j . 
it was assumed that all neutron stars have the same radii 
and they should lie between 7.6 and 10.4 km and 10.9 and 
12.7 km respectively. We resume the properties of 1.4M0 
quark stars obtained with the formalisms discussed in the 
present work in Table|2| We see that for all magnetic fields, 
the MIT bag model fulfills all constraints, except the ones 
proposed in refs. [5S] and [55] . because they are mutually 
exclusive and cannot be satisfied at the same time. 


3 The NJL model 

The EOS of magnetized matter obtained with the NJL 
model has already been extensively discussed in [^ , [23] , 
[55] , [57] . The Lagrangian density is the same as given for 
the MIT bag model in Eq. m and the leptonic sector is 
also the same as described in Section[2] The quark sector is 
described by the su(3) version of the Nambu-Jona-Lasinio 
model 

Cf = Ipf [7^ {id^ -qfA^)- Ac] V'/ + Bsym + Bdet ■ 

( 20 ) 
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Table 2. Properties of the IAMq quark star 


Model 

B 

^^max 

^^max 

R 

Sc 


(G) 

(Me) 

(Mq) 

(km) 

{im-Y 

EEIT isotropic 

TcF 

1.40 

1.71 

10.23 

2.64 

EEITanisotropic 

TcF 

1.40 

1.72 

10.24 

2.54 

MITchaotic 

TcF 

1.40 

1.72 

10.24 

2.65 


EEIT isotropic 

10 "“ 

1.40 

1.72 

10.25 

2.65 

EEITanisotropic 

TcP 

1.40 

1.72 

10.24 

2.65 

MITchaotic 

TcP 

1.40 

1.73 

10.25 

2.65 


EEITisotropic 

3 . 10 "“ 

1.40 

1.73 

10.34 

2.53 

EEIT anisotropic 

3 . 10 "“ 

1.40 

1.73 

10.13 

2.78 

MITchaotic 

3 . 10 "“ 

1.40 

1.74 

10.32 

2.60 


EEITisotropic 

10""' 

1.40 

1.81 

11.10 

2.07 

EEIT anisotropic 

TcP 

1.40 

1.78 

9.85 

2.97 

MITchaotic 

TcP 

1.40 

1.81 

11.09 

2.09 


The Lsym and Cdet terms are given by: 

8 

^sym — ^ [W/AaV'/)" + (^/*75Aa^/)"] , (21) 

a—0 

Cdet = -K {det/ + det/ [■0/(1 - 75 )'*/'/] } 

( 22 ) 

where ipf = (tt,d, s)^ represents a quark field with three 
flavors, ihc = diagj(mu, nis) with niu = nid yf mg 
is the corresponding (current) mass matrix while qj = 
diag((7u, Qd, 9s) is the matrix that represents the quark 
electric charges. Aq = with I being the unit ma¬ 

trix in the three flavor space, and 0 < Aa < 8 denote 
the Gell-Mann matrices. The term {Cdet) is the t’Hooft 
six-point interaction and Cgym is a four-point interaction 
in flavor space. The model is non renormalizable, and as 
a regularization scheme for the divergent ultraviolet inte¬ 
grals we use a sharp cut-off A in three-momentum space. 
In the present work, we use the HK parametrization pro¬ 
posed in [B3] : A = 631.4 MeV , m„ = rud = 5.5 MeV, 
rus = 135.7 MeV, GA^ = 1.835 and KA^ = 9.29. 

The conditions of charge neutrality and /3-equilibrium 
given in Eqs. (nni) and m are also enforced. 


3.1 NIL - Isotropic EOS 

In the mean field approximation the pressure can be writ¬ 
ten as 

Pf = &u+Sd + 0s — ‘^G{4i‘^ + (j)d+4’'^)+'iK(j)u4‘d4‘s I (23) 

where the effective quark masses can be obtained self con¬ 
sistently from 

Mi = rrii - 4G(j)i+ 2K4>j4>k, (24) 


with {i,j,k) being any permutation of {u,d,s), 


9f = ( 0 ^° -b -b ^ , (25) 


where the vacuum contribution reads 




{A + Cyl) 
Wf 


-€aA (A^ + e^) I , 


(26) 

and with ca = \J + M'j , the finite magnetic contribu¬ 
tion is given by 


1 

C'(-l,a^/) - i^{x) -Xf)\nxf + -X 

(27) 

with a;/ = M'j / {2\qf\B) andC^(—l,x/) = dC,{z,x j) / dz\z=- 
where C,{z,Xf) is the Riemann-Hurwitz zeta function. The 
medium contribution can be written as 


gjag ^ 


Nc{\qf\Bf 

27r^ 


nmed 


E 


O^k 


\qf\BNe 

Att'^ 


Sf{k,BY^ 


—Sf{k, B)^ In 



- Sfik,B)^ 
Sf{k,B) 


(28) 


where Sf{k,B) = y + 2|9/|ilfc, and the upper Landau 
, level (or the nearest integer) is defined by 


fil — Mf 

2\qf\B 


(29) 


The condensates (/)y in the presence of an external mag¬ 
netic field can be written as 


where 


= - 


= + + (30) 

A + eA 


MfNe 

27r^ 




AeA — Mj In 

M}\qf\BNe 

2'k‘^ 


Mf 


, (31) 


lnT(a;/) - ^ ln(27r) +Xf-^ {2xf - 1) ln(x/) 


(32) 


and 


—j 

(33) 

The corresponding leptonic contributions are the same 
as the ones used for the MIT bag model. The final expres¬ 
sions for the pressure and energy density are given as in 
Eq. ([T^ and p3|l . 


^jed 
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3.2 NJL - Anisotropic EOS 

The parallel and the perpendicular components of the 
pressure can be written in terms of the magnetization as 
in Eq. (HI- The complete calculation is shown in m. but 
the main formulae are given next. 

For the quark sector the derivatives of the pressure 
with respect to the magnetic field are: 


dPt 

Mf = “ ^G{(l)u4’u + + 4’s(l}'s) 


B = 1x10 
B = 1x10 °G 
B = 3x10 “G 
B = 1x10^®G 


1 - 1 - 1 -r 


cs 


W) 

c3 


where 

and 

with 


I- AK{4>'^4>d4>s + 4‘u(l)'d4>s + 4>u(t>d(t>'s) 

r\/ (nt vac , nf mag . n/med\ 

Uf = [Uj: +e/y )Mf , 

0/ _ vac I ‘^CLg I med\ 



^,pag 


>/ 


)Mf , 


(34) 

(35) 

(36) 


0.2 0.4 0.6 0.8 


e‘ 


/ mag _ 2 


NMfWj 


B 


47r^ 


lnr(a;/) - ^ln(2T:)+Xf - [xf - i)ln(a:/) 


nmed at nl I ^ max 


(37) 


B 


27r^ 


i/mag 


fe=0 



N^Mj 


B 


Stt^B 

{\Qf\B + Mf iP<^°\xf) -\n{xf) I , 


(39) 


where ip^i^Xf) = ^rlxp digamma function. The in 

medium contribution reads 


'led 


if med 
<?/ = 


B 


- + 


Nc\qf\B 

27r2 


E 


Oik- 


/r/M/(fc|g/| 


=0 Sf{k,B)^Jn‘j-Sf{k,B) 


(40) 


and and vanish. 

For the leptonic sector one easily gets 

^pmed pmed 

Ml = ' — * 


dB 


B 




ik\m\) , (41) 


fc=0 


which is another way of expressing Eq. (1151) with the ap¬ 
propriate substitutions for the leptons. 

The total parallel and perpendicular pressures are given 
as in Eq. (USD- 


1 1.2 1.4 1.6 1.8 

n (fm'^) 


Fig. 12. Magnetization as a function of the number density 
for the NJL model. 


3.3 NJL - Chaotic field 

Again, as in the MIT model, the final expressions for the 
EOS are given by Eas. dT^ and (ITOl) . 


3.4 Results NJL 


We start by showing the EOS obtained from the three for¬ 
malisms within the NJL model. We have used eo = 7.81 
fm“^ in Eq. © because this value is the central energy 
density of the maximum mass obtained with the NJL 
model for non-magnetized matter with the HK parameter 
set [03] ■ As in the MIT case, the magnetic field is taken 
as constant in Pf and Pi and Eq.(|T|) is used in the terms 
proportional to B^ only and the same holds for the energy 
density terms. 

In Fig. ll2l we show the magnetization of the system for 
different values of the magnetic field. As already discussed 
in [67] and also seen in Fig. [TJ the number of van alphen 
oscillations, related to the number of spikes in the NJL 
model, are larger for lower magnetic fields and decrease 
for stronger magnetic fields, when the number of filled 
Landau levels is small. Whenever fj,f = s/, i.e., whenever 
B approaches a n 7 ^ 0 Landau level, the denominator in 
Eq. (H01) becomes zero and the contribution of this term to 
the magnetization of the system generates the spikes. This 
kind of contribution is not present in the MIT model and 
hence, the overall pattern is quite different in both models. 
Another important difference refers to the lack of oscilla¬ 
tions in between 0.5 and 0.7 fm“^ for all the fields within 
the NJL model. As pointed out in m, this behavior is due 
to the u and d quark restoration of chiral symmetry and 
before the onset of the s quark. Again, this behavior is not 
seen in the MIT model, where chiral symmetry effects are 
not present. Finally, it is worth mentioning that the MIT 
bag model only produces negative magnetization for very 
low magnetic fields while the NJL model gives negative 
values for the magnetization for all fields. Bearing in mind 
that negative magnetization refers to magnetic repulsion 
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Ph 


2.4 
2 

1.6 
1.2 
0.8 
0.4 
0 

012345678 

8 (fm'^) 


^ I I I I I I 

P(perp) - 

P{parall) 

B = IxlO^^G 


I l.ll .A><' 



Fig. 13. NJL model EOS obtained for B = 10^^ G without Fig. 14. NJL model EOS obtained for B = G without 
the term proportional to B^ the term proportional to B^ 


(diamagnetism) as opposed to positive magnetization that 
referes to magnetic attraction (paramagentism), the MIT 
and the NJL models show different physical pictures. We 
address this point again in the next section. 

In Figs. [HI [Til [T51 and m we show the anisotropic 
EOS with parallel and prependicular pressures without 
the inclusion of the term and in Figs. [T71 [THl [12] and 
HDl the same EOS are plotted with the inclusion of this 
term. Eor B — 10^^ G, the contribution of the B^ term is 
almost negligible and hence, we can see that the magneti¬ 
zation already plays a role in the pressures, which are no 
longer coincident, as in the MIT model. As B increases, 
the effects of the magnetization are clearly noticed. Once 
more, as in the MIT model, if the B^ contribution is taken 
into account, the parallel pressure increases at lower en¬ 
ergy densities and then decreases and goes to zero, reach¬ 
ing this point at lower values for larger magnetic fields. 
However, if one compares the isotropic pressure with the 
perpendicular pressure, they are just similar in average, 
since the magnetization effects are clearly present. This 
behavior is also quite different from the one exhibited by 
the MIT model and discussed in a previous section. Pre¬ 
cisely because of the discontinuities in the perpendicular 
pressure, the anisotropic EOS cannot be used as input to 
the TOV equations. 

At this point it is worth mentioning that similar calcu¬ 
lations were perfomed in [68] for the NJL model with and 
without a vector interaction. However, in [68] the spikes 
in the EOS are not seem for fields as large as 4 x 10^® G 
. We have checked that if we plotted the pressure and en¬ 
ergy density in MeV/fm® (without the B^ term) as in the 
mentioned reference, the spikes are almost imperceptible 
in the region of energy densities corresponding to 500-1000 
MeV/fm®, but they are indeed present since they appear 
whenever /rj = s/ in the denominator of Eg. (1401) . Notice 
also, that the parametrizaton of the electromagnetic in- 
teration in the EOS (exactly the B^ term) is different in 
both papers, resulting in different quantitative results. If 
this term grows too rapidly the oscillation-free zone can 
dominate the equation of state. 



Fig. 15. NJL model EOS obtained for H = 3 x 10^® G without 
the term proportional to B^ 



Fig. 16. NJL model EOS obtained for B = 10^® G without 
the term proportional to 
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Fig. 17. NJL model EOS obtained for B — 10^^ G with the 
term proportional to 



Fig. 18. NJL model EOS obtained for B — 10*^® G with the 
term proportional to B^ 



Fig. 19. NJL model EOS obtained for B = 3 x 10*^® G with 
the term proportional to B^ 



Fig. 20. NJL model EOS obtained for B — 10^® G with the 
term proportional to B^ 



Fig. 21. NJL EOS for different values of the chaotic magnetic 
field. 



Fig. 22. NJL model: Mass-radius relation within isotropic 
pressure for different values of magnetic field. 
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Fig. 23. NJL model: Mass-radius relation within chaotic mag¬ 
netic field approximation for different values of magnetic field. 


Table 4. Properties of the IAMq quark star 
Model B Mmax Mb max R £c 

(G) (Mo) (Mo) (km) (fm-4) 

NJLisotropic icp Oo L46 TOS IM 

NJLchaotic ICF L40 L46 TOS IM 

NJLisotropic icP Oo 06 90 Ira 

NJLchaotic IcP Oo 06 T05 l82 

NJLiaotropic IT(P Oo 06 T05 l60 

NJLchaotic 3.10"“ Oo 06 T04 ISO 

NJLiaotropic IcP Oo Os 8(9l los 

NJLchaotic IcP Oo 03 Ptt OI 


Table 3. Properties of the maximum mass quark star 


Model 

B 

^'^max 


R 



(G) 

(Me) 

{Me) 

(km) 

(fm-h 

N JLisotropic 

10 "'■ 

1.45 

1.53 

8.93 

7.51 

NJLchaotic 

10 "" 

1.45 

1.53 

8.92 

7.66 


N JLisotropic 

10 "“ 

1.46 

1.54 

8.87 

7.76 

NJLchaotic 

10 "“ 

1.45 

1.53 

8.92 

7.46 


NJL 

isotropic 3.10"“ LSS L62 Oi 058 


NJLchaotic 3.10"“ Ot L54 085 7.86 

Ml 

isotropic OcP Oo L77 02 10.25 

NJLchaotic IcP 09 09 OI 8.93 



n (fm"^) 


Fig. 24. Magnetic field times the system magnetization for 
B = 3 X lO"® G for the MIT and the NJL models. 


Finally, in Fig. [511 the EOS obtained with the chaotic 
formalism is shown for different values of the magnetic 
field. As the magnetic field increases the EOS becomes 
harder and within the NJL model, they are more sensible 
to the magnetic field than within the MIT model, where 
the EOSs just differ from each other for fields of the order 
of 10^® G, as seen in Fig. [8] 

We then use the EOS obtained with the isotropic and 
chaotic formalisms to compute the stellar properties and 
the results are displayed in Figs [22l and [23] respectively, 
and in Table [3] we summarize the properties of the maxi¬ 
mum mass quark star. As already expected, from the re¬ 
sults existing in the literature [53], [30], the maximum 
stellar masses are very low, indicating that the NJL model 
cannot explain very massive magnetars. A solution to this 
setback is to include the vector interaction in the NJL 
model as in [30] . The overall conclusions coming out of the 
macroscopic properties for the NJL model are the same as 
the ones obtained with the MIT model, i.e., the maximum 
mass increases much more from a non-magnetized star to 


a strongly magnetized one if the isotropic formalism is 
used than if the chaotic field approximation is assumed. 

Finally we display the properties of the canonical 1.4M0 
neutron stars in Table S) We see that with the NJL, the 
quark star radii are very low, and are always in accordance 
with ref. [BS] , instead the constraints of ref. [BB] . 


4 MIT versus NJL - differences 

We start by reanalyzing the effects of the magnetization 
in the perpendicular pressure of both models. According 
to Eq. (II6L the magnetization always comes multiplied 
by the magnetic field. We then plot, in Fig. |M| the term 
that really influences the perpendicular pressure for B = 
3 X 10^® G. Theses results help us to understand the effects 
observed in Figs. |3| and [T3| for the MIT and NJL models 
respectively. The term proportional to the magnetization 
is very small in both models and even when it becomes 
negative, the perpendicular pressure is not much smaller 
than the parallel one. 
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B (G) 


Fig. 25. Anisotropic pressures for the MIT and the NJL with¬ 
out the contribution for n = 0 .2no. 



Fig. 26. Anisotropic pressures for the MIT and the NJL with¬ 
out the contribution for n = 5no- 


Our last figures, Fig. [^S] and Fig. [HI show the perpen¬ 
dicular and parallel pressures obtained from both models 
at a fixed densities, 0.2no and 5no with uq = 0.16 fm“^ be¬ 
ing the nuclear matter saturation density. These densities 
are chosen to represent stellar matter close to the crust 
and in the core of the star respectively. In these graphs 
the terms proportional do are not included. As shown 
in [41] for the MIT model, both pressures do not differ 
much up to a magnetic field of the order of 10^® G, when 
they start to deviate. For larger magnetic fields, of inter¬ 
est, for instance, in heavy ion collisions, the deviation first 
increases and then decreases again, stabilizing around 10^° 
G. For the NJL model, the behavior is not very different 
in average. For the case of lower density, the perpendicular 
pressure presents the typical spikes generated by the term 
that includes the magnetization and the parallel pressure 
is lower than the corresponding parallel pressure obtained 
with the MIT model, but the general trend is the same 
in both models. For n = 5no, the number of spikes ob¬ 
tained with the NJL model decrease because this density 
corresponds to the region where the magnetization is very 


low, as shown in Fig. [T^ One can see that the larger the 
density, the larger the differences between the pressures 
for fields higher than I0^®G. 

Another topic of interest is the effect of the magnetic 
field on the quark star radii. We can also see from Tables 2 
and 4, that the NJL model always yields lower maximum 
masses and lower radii for the canonical IAMq NS. Com¬ 
paring Figs, ini and |TT] with Figs. |H1 and |H1 we see that 
while the magnetic field increases the maximum mass in 
both MIT and NJL models, the effect on the radii is op¬ 
posed. In both formalisms, isotropic and with chaotic field, 
for a fixed mass, the radii increase with the increase of the 
magnetic field in the MIT bag model and decrease in the 
NJL. In fact, the only case that produces a more com¬ 
pact quark star in the MIT bag model is when anisotropic 
pressure is considered. Unfortunately we are unable to cal¬ 
culate the TOV equation with the anisotopic NJL EOS 
due to the high oscillations in the pressure. Nevertheless, 
even if we could, these results would have to be taken with 
care since, as discussed in Section [5] because just one of 
the pressures is taken into account, since the other one 
becomes zero at still quite low densities. To trust macro¬ 
scopic results obtained from a TOV-like equation, a self- 
consistent model for the study of the quark star structure 
such as the one proposed in [TO] should be done, so that 
more precise conclusions could be drawn. 

5 Final remarks 

In the present work we have used two models widely used 
to describe quark matter, the MIT and the NJL model and 
checked the effects of strong magnetic fields on their mag¬ 
netization, EOS and stellar properties. To accomplish this 
task, we have revisited three formalisms generally applied 
to describe magnetized matter: the first one assumes that 
the EOS is isotropic, the second one that it is anisotropic 
and the third one that a chaotic field is possibly generated 
and maintained in the system. 

For the second formalism, the calculation of the mag¬ 
netization of the system is mandatory, since one of the 
pressures depends on this quantity. However, independently 
of the formalism used, the magnetization is a quantity that 
should be understood in a magnetized system. We have 
verified that both models produce quite different results. 
The NJL model presents three terms for the magnetiza¬ 
tion coming from the medium, the vacuum and the mag¬ 
netic field itself. Just the first one is present in the MIT 
model and this explains most of the differences between 
both models. 

As for the EOS, the isotropic and the chaotic for¬ 
malisms present similar patterns in both models, the NJL 
being always slightly more sensible to the magnetic field 
strength. However, the anisotropic EOS shows a quite dis¬ 
continuous behavior due to the spikes that appear in the 
calculation of the magnetization of the system. 

When the macroscopic properties are computed, the 
isotropic EOSs always provide maximum stellar masses 
that are much larger than the ones obtained with a non- 
magnetized matter, in contrast with the values obtained 
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with the chaotic field that do not result in much higher 
maximum masses. The small increase in the maximum 
masses is also found in more sophisticated calculations 

EH, IZO]. 

Other important aspects related to quark matter sub¬ 
ject to strong magnetic fields are the polarization and the 
viscosity of the system. For totally polarized matter, all 
particles lie on the lowest Landau level. The quark spin 
polarization has already been studied for both the MIT 
bag model in m and the NJL model in im, m and 
hence, we only comment on the already known results. 
The polarization of the system increases with the increase 
of the magnetic field. In m, it is shown that the total 
polarization obtained with the MIT bag model occurs for 
i? ~ 2 X 10^® G for matter in /3-equilibrium. This result 
depends on the density and the temperature considered 
and is different for each quark flavor, since it depends on 
the charge of the quark m- For matter in /3-equilibrium 
described by the NJL model, the electron becomes totally 
polarized for magnetic fields larger than 9 x 10^^ G and the 
s quark before 10^® G. The u and d quarks become totally 
polarized for approximately the same magnetic field as in 
the MIT bag model, i.e., higher than 10^® G and lower 
than 10®° G, the exact value being parameter dependent. 

For quark matter, the non-leptonic processes 

u + s u + d, u + d u + s 

are responsible for the restoration of chemical equilibrium 
and hence, they determine the bulk viscosity of the sys¬ 
tem. In m, bulk (and shear) viscosities are derived and 
discussed. Similar calculations remain to be done with the 
MIT bag model for the chaotic field approximation and for 
the NJL model and the results can shed some light on the 
stability of magnetized quark matter. 

The physics of quark star radii is another open puzzle. 
It is well known that the radii of hadronic neutron stars 
are correlated with the symmetry energy slope L, |72] . 
m, m, m, m- However we do not know yet if there 
is an analogue physical parameter that controls the quark 
star radius. 

As mentioned in the Introduction, the influence of the 
quark AMM should not be completely disregarded. In 
it was shown that for quark matter the scale for the pertur¬ 
bative approach is defined by the constituent quark mass 
and therefore, the effects of the AMM should also be con¬ 
sidered in the description of magnetized quark matter. As 
the quark mass is about one third the nucleon mass, the 
critical field will be approximately one order of magnitude 
smaller than the nucleon critical field, but still larger than 
the maximum field expected inside a quark star. Hence, 
the degeneracy of Landau levels will certainly be affected 
by the inclusion of the coupling between the field and the 
fermion anomalous magnetic moment. In a future work 
this problem should also be considered. 

Before we conclude, it is important to stress that we 
have not exhausted the discussion on all formalisms since, 
as already mentioned in the Introduction, there are still 
other ones in the literature. One possible approach to 
deal with the anisotropic pressures and avoid the use of 


the TOV equations, is to follow the prescription used in 
[31] with the help of the Hartle-Thorn method. We have 
seen that effects of the magnetization should not be disre¬ 
garded as far as magnetic fields are stronger than 10^^ G 
and hence, the calculations performed in EH, which com¬ 
pletely ignore the contribution due to the magnetization of 
the system in the perpendicular pressure, should be revis¬ 
ited. There are also more sophisticated calculations that 
consider the anisotropy in solving Einstein’s field equa¬ 
tions in a fully general relativistic formalism [55] , [56] , [70] , 
but the computational price paid is huge. Moreover, the 
fields generated at the surface of the star in EO] are of the 
order of 10^^ G, much higher than the 10^® G magnetic 
field expected in the surface of magnetars. 
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